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Motivated by recent experiments [A. Lorke et al., Phys. 
Rev. Lett. 84, 2223 (2000)] an analysis of the ground state and 
far-infrared absorption of two electrons confined in a quantum 
ring is presented. The height of the repulsive central barrier in 
the confining potential is shown to influence in an important 
way the ring properties. The experiments are best explained 
assuming the presence of both high- and low-barrier quantum 
rings in the sample. The formation of a two-electron Wigner 
molecule in certain cases is discussed. 

PACS 73.20.Dx, 72.15.Rn 



I. INTRODUCTION 

The possibility of confining a controlled number of 
electrons in artificial nanostructures at the interface of 
two semiconductors, achieved in recent years ||^, has at- 
tracted much interest. An obvious reason is the poten- 
tial application of these systems to sub-micron electron- 
ics technology but, from the theoretical point of view, 
much importance is attributed to their character of novel 
quantum systems in which new physics, as well as analo- 
gies with atomic clusters, atoms or nuclei, may be found. 
The most widely studied systems have been the so-called 
'quantum dots', in which a certain number of electrons 
form a two-dimensional compact island, modelled by a 
confinig potential of parabolic type (for small dots) plus 
the mutual electron-electron interactions. The electronic 
islands have been also formed in the shape of rings. In 
fact, mesoscopic rings were produced by Dahl et al. and 
their magnetoplasmon resonances measured for different 
ring widths Subsequent theoretical models, using 
semiclassical methods |^ and density functional theory 
Q provided a good interpretation of the measured ring 
spectra. There also exist approaches based on exact diag- 
onalization methods for very small numbers of electrons 
which have addressed the so called 'persistent current' 
1^ and the optical excitations for very narrow rings (ap- 
proaching the Id limit) as well as for dots with a repulsive 
scatterer center 1^,0. Besides, a description of the opti- 
cal properties in two electron rings based on an analytical 
treatment of the ground state was given by Wendler et 
al. §. 

Very recently, Lorke and co-workers reported the fab- 
rication, using self-assembly techniques, of electron rings 
in the true quantum limit, not affected by random scat- 
terers, containig two electrons M. The far-infrared (FIR) 



spectra as a function of the magnetic field has shown the 
appearance of electronic excitations peculiar of a ring ge- 
ometry as well as a magnetic-field induced transition in 
the ground state angular momentum. The two electron 
system is the smallest one with electron-electron interac- 
tion and its reduced size places it within the reach of es- 
sentially exact methods, such as the Hamiltonian diago- 
nalization in a basis of functions mentioned above. In this 
paper we have followed an alternative although also exact 
method, based on the solution of the Schrodingcr equa- 
tion by discretizing the coordinate space in a uniform 
grid of points. The method is therefore free from basis 
selection and truncation although its limitation lies, ob- 
viously, in the number of spatial points. This scheme can 
also be used to obtain the dynamical properties by solv- 
ing the time-dependent Schrodinger equation and moni- 
toring the evolution of relevant physical quantities with 
time. In this work we have addressed the linear response 
regime, for small amplitude oscillations. However, one of 
the potentialities of the real-time approach is the study 
of non- linear (large amplitude) dynamical processes. By 
finding the exact solution one can easily quantify the ef- 
fect of the electron-electron interaction as well as check 
the validity of approximate descriptions, such as density- 
functional theory. In fact, an analysis of the experiment 
based on the local spin-density approximation (LSDA) 
in a symmetry-restricted approach has already been pre- 
sented in Ref. ]lO[ |. 

An important input to the theory is the external con- 
fining potential acting on the electrons. Since we aim at 
a direct comparison with experiment, the strategy should 
be to introduce a reasonable model potential, depending 
only on a few parameters and explore the exact solution 
as a function of this parameters. Our approach has been 
to take the parabola parameter ujq and ring radius Rq as 
guessed in Ref. |9| and explore the effect of the central 
repulsive barrier by solving for a low-, an intermediate- 
and a high-barrier ring. With these barriers the system 
changes from a compact dot-like structure to a well devel- 
oped ring one. We show that the barrier height modifies 
the ring properties and that the experiments are best ex- 
plained assuming a mixed sample containing both low- 
and high-barrier rings. 

In recent papers Yannouleas and Landman , as well 
as Koskinen et al. |12| have discussed the electronic lo- 
calization, in the relative frame, known as the forma- 
tion of Wigner molecules in two electron quantum dots 
and in rings with 6 electrons, respectively. Both groups 
have addressed the problem by calculating the low-lying 
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states and discussing their rotor-like structure. This crys- 
taUization is in good agreement with the prediction of 
mean-field (MF) theories like the Hartree-Fock or the 
Local-Density-Functional approaches. We also comment 
on the formation of a Wigner molecule in the high-barrier 
rings by comparing the conditional probability distribu- 
tion with the LSDA density. In this case, the LSDA pre- 
dicts a symmetry-broken ground state that nicely agrees 
with the electronic distribution built from the condi- 
tional probability density of the exact solution to the 
Schrodinger equation. 

The structure of the paper is as follows. Section II is 
devoted to the exact solution for the ground state and 
Sec. Ill to the corresponding exact solution for the dy- 
namics. In Sec. IV we discuss the comparison with the 
mean field approach and the formation of the Wigner 
molecule. Finally, the conclusions are presented in Sec. 
V. 



II. GROUND STATE 

A. Hamiltonian and resolution method 

The Hamiltonian of the two-electron system, in a mag- 
netic field, using effective units is 



i(-iV,+7A(r,))'-|-y(r,) 



i=l,2 

^ +.g*m*^B.S 
ri-r2 2 



(1) 



where A is the vector potential, 7 = e/c (we assume 
Gaussian magnetic fields) and V{t) is the confining ex- 
ternal potential. The electronic positions are restricted 
to the xy plane, i.e., r = {x,y), and we consider a con- 
stant and perpendicular magnetic field B = Be^, de- 
scribed in the symmetric gauge with a vector potential 
A(r) = B/2{—y,x). The last piece is the Zeeman term, 
depending on the total spin S = Si + S2 the effective 
gyromagnetic factor g* and electronic effective mass m* 

0- . 

The ring external potential is modelled by the following 
circularly symmetric piecewise function. 
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if r > i?o 



(2) 



c(r - i?o) + d{r - Roy] if r < i?o 



In this expression Rq is the ring radius and luq gives the 
curvature around the minimum at Rq. Fixing a barrier 
height Vfc at the origin determines the parameters c and 
d, with the additional requirement of smootheness at Rq 

Since Hamiltonian (pi) depends on spin only through 
the Zeeman term, the spin part of the wave function 



is given by the well known singlet and triplet combina- 
tion of spinors. Our method to obtain the spatial wave 
function is based on the only assumption that its gra- 
dient and Laplacian can be obtained by using finite dif- 
ference formulas in a spatial grid of uniformly spaced 
points. Associating a macroindex / = (xi, yi, X2, j/2) to 
each point and dicretizing the derivative operators (we 
use 5 point formulas), the Schrodinger eigenvalue equa- 
tion transforms into the linear problem 



(3) 



It is worth to mention that the interaction and poten- 
tial terms are diagonal (local) and only the kinetic con- 
tribution gives non diagonal elements. Besides, since the 
derivatives involve only a finite number of points (we have 
tipically used 5-point formulas) the Hamiltonian matrix 
is very sparse. We have solved Eq. (^ for the ground 
state using the iterative imaginary time-step method [p^ . 
In principle, one can also obtain excited states by intro- 
ducing the constraint of orthogonalization to the lower 
eigenvectors, although the computational effort increases 
very rapidly. A very useful optimization is introduced 
by considering only symmetric (singlet) or antisymmet- 
ric (triplet) wave functions with respect to the exchange 
of the two particle positions, since this allows to reduce 
by a factor two the wave function dimension. 

The above method has proved to be quite stable and 
the convergence with the number of points very fast. 
Typically, we have used from Ni «20 up to «50 points 
for the discretization of a single dimension. Notice that 
the wave function dimension (iVj*) would normally im- 
ply that the Hamiltonian matrix {N* x N*) exceed the 
storage capability. This is avoided by fully exploiting the 
sparseness of the Hamiltonian and not really dimension- 
ing such a big matrix. Analogously, the imaginary-time 
step method allows us to determine the lower eigenstate 
of this matrix which would be unfeasible by direct di- 
agonalization methods. We have also checked the con- 
vergence with the number of points used to discretize 
the Laplacian operator, e.g., when going from 5 to 9 
point formulas a relative difference less than 1/10'^ in 
the ground state energy is found. 



B. Results 

In this subsection we present the results for the ground 
state of rings having Rq = 1.4 Oq and ojq = 1.1 H*. The 
corresponding physical values are, approximately, 14 nm 
and 13 meV and are taken as reasonable values from the 
experiment. Three different barrier heights are consid- 
ered Vb = 0.75, 3 and 12 H*. As shown in Fig. 1 these 
three values correspond, going from low to high barrier, 
to decreasing central densities. In what follows the three 
rings will be referred to as VI, V2 and V3 for increasing 
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barrier height, respectively. While ring VI has a compact 
density because of its rather low barrier, V3 has a real 
ring shape, with a hollow central region. 

Figure 2 displays the expectation value of the angular 
momentum (L^) as well as the energy of the lowest level 
for both singlet and triplet states as a function of the 
magnetic field. At i? = the ground state is a singlet 
and with increasing B there is an alternation of triplet 
and singlet ground states in the three rings. The first 
singlet-triplet transition occurs approximately at 3, 2 and 
1 T for VI, V2 and V3, respectively. The triplet-singlet 
energy differences remain quite small after the first cross- 
ing which hints at possible transitions induced by small 
non-spin conserving perturbations like spin-orbit inter- 
actions or magnetic impurities. We also notice how the 
ground state increases in steps its angular momentum (in 
absolute value) as a function of B. This reflects the ex- 
act property that angular momentum is a good quantum 
number since the Hamiltonian commutes with the op- 
erator. The angular momentum gain with magnetic field 
is faster for the higher barrier rings. In fact, for V3 the 
evolution is markedly linear, already indicating a rotor- 
like evolution for this ring. This is further discussed in 
Sec. IV. 



C. Non-interacting solution 

When neglecting the Coulomb interaction the Hamil- 
tonian (|l|) separates for each particle and, therefore, a 
solution built from single-particle (sp) orbitals with ra- 
dial n and angular £ quantum numbers is adequate [Q. 
Figure 3 shows the sp eigenvalues for all, low and high 
barrier, rings. In fact, from the lowest maximum of 
the capacitance-voltage spectrum Lorke et al. obtained 
the single-electron ground-state shift with magnetic field. 
This manifests a change in slope at -B ^ 8T that was ex- 
plained as the transition from £ = 0to^=— lina model 
potential. Of course, this is in agreement with the sp lev- 
els for VI, since the low-barrier ring essentially coincides 
with the model potential used in Ref . . The first cross- 
ing for V2 is located around 4T, while as shown in Fig. 3 
the high barrier case (V3) has again a change in slope at 
B ~ 8T, attributed to the £ = -1 — > £ = -2 transition, 
which may indicate that the experimental sample could 
contain both low and high barrier rings. This point will 
be further discussed when presenting the FIR spectra. 

The non-interacting ground state with two electrons is 
obtained by filling the lowest orbitals of the sp scheme, 
with the constraint of having a singlet or triplet spin. 
The corresponding plots for the evolution of (L^) and 
the ground state energy as a function of magnetic field 
are displayed in Fig. 4. By comparing it with Fig. 2 we 
notice the effects of the Coulomb interaction between the 
two electrons. Obviously, the non-interacting case gives 
lower ground state energies since the Coulomb repulsion 



is neglected. Due to the increasing effect of the con- 
fining potential the difference is higher for VI than for 
V3. Differences in angular momentum are also higher 
for the low barrier ring. This is indicating that in well 
developed rings the interaction effects are less important 
than in more compact structures, in agreement with the 
findings of Ref. Q and with the obvious expectation of 
a larger VJ, keeping the electrons more apart from each 
other than a lower barrier. 

We also realize that the triplet solution is unfavoured 
in the non-interacting case. In VI the gs has always 
5 = 0, except for i? = 8 T for which singlet and triplet 
are essentially degenerate. A similar thing happens in the 
other rings, with the difference that the singlet-triplet de- 
generacy points appear at lower B's. These differences 
with respect to the exact results, in which clearer singlet- 
triplet oscillations are seen, indicate the importance of 
the Coulomb interaction for a proper ground state de- 
scription. 



III. TIME INTEGRATION 

The time evolution of the two-electron wave function 
is given by the time-dependent Schrodinger equation 



l—^{Yi,Y2\t) 



(4) 



where H is the Hamiltonian operator (|^). We have solved 
Eq. (Q), in the same spirit as for the ground state, by 
discretizing the coordinate space in a uniform grid of 
points. Time is also discretized and the evolution from 
time (n) to {n+1) is obtained through the unitary Crank- 
Nicholson algorithm 
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This is an implicit scheme, in which the new v]/("+i) may 
be obtained by iteration. The above algorithm allows 
to perform the dynamical simulation of Eq. (^) by using 
a small time step At. Of course, the feasibility of this 
approach to the solution of Eq. relies on the very small 
number of constituents (two electrons) of the system. 

In this paper we discuss the FIR absorption of two- 
electron rings and therefore our interest will focuss on the 
dipole states corresponding to small amplitude charge os- 
cillations. The ground state wave function is perturbed 
by an initial (small) rigid displacement in a certain direc- 
tion e, representing the electric field direction, and the 
dipole moment evolution is then followed in time. The 
absorption cross section is given by the frequency trans- 
form of the dipole signal 
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dt(^'(i)|£'|*(t)) exp{iojt) 



(6) 
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where the dipole operator is D = e ■ (ri + r2). As a 
test of the method, developed from previous calculations 
within mean-field approaches, we have checked that the 
oscillation frequencies when assuming a pure parabolic 
potential are given, as a consequence of the generalized 
Kohn theorem, by the values 

uj± = ^ujI+lj^JA ± Ljj2 , (7) 
where lUc = eB/c is the cyclotron frequency. 

A. FIR spectra 

Figure 5 shows the absorption spectra for singlet and 
triplet states corresponding to different magnetic fields. 
The intersection of the baseline for each spectrum with 
the horizontal axis indicates the magnetic field in each 
case. Notice that the horizontal scale for each spectrum 
is arbitrary, and therefore, it is not given by the hori- 
zontal axis. The spectra are drawn in linear scale and 
arbitrary units while the vertical scale gives the energies 
in meV. Not surprisingly the VI spectrum is very similar 
to that of compact dots, with two branches of opposite B- 
dispersion which merge at i? = 0. This is indeed evident 
by comparing the spectra with the analytical expressions 
(dashed lines), where an empirical luq value fixed to 
the lowest B = peak energy has been used. Except for 
some fragmentation details these two branches are rather 
similar for both singlet and triplet cases. More impor- 
tantly, we want to stress that essentially no strength is 
contained in the energy region 15 meV < a; < 25 meV 
for S < 6 T. This fact does not imply that there are no 
excitations in this region. In fact they can be seen using 
logarithmic scale, as has been done in the LSDA approach 
of Ref . . We have sticked however to the linear scale 
firstly because the real-time approach requires extremely 
long simulation times in order to extract so low intense 
frequencies, and secondly we do not think that the ex- 
periments can easily identify such low intensity peaks. 

The V2 spectrum and, more clearly, the V3 one reflect 
a behaviour qualitatively different from that of VI. The 
high energy branch does not merge with the low one for 
vanishing B, deviating from the analytical prediction of 
Eq. . This behaviour is well known for quantum rings 
and is a manifestation of the violation of Kohn's theo- 
rem in these systems [|[j8). For V2 a global gap w 8 — 14 
meV persists at all magnetic fields for both spins, while in 
the case of V3 the corresponding gap lies approximately 
in the region 7—17 meV. Therefore, this gap consti- 
tutes a direct manifestation of the barrier height in the 
FIR response. The B-evolution of the spectra is char- 
acterized by small jumps in the peak positions for each 
branch which must be attributed to changes in total an- 
gular momentum. Since V3 is the one gaining more an- 
gular momentum its spectrum also displays more abun- 



dant discontinuities. We also notice that when the ring 
barrier grows (from V2 to V3) an important amount of 
strength is placed in the interval 15-25 meV at low mag- 
netic fields. Actually the V3 peak in this region has a 
strength comparable to the one at low energy w 50%, as 
found experimentally [|lO| , while our calculation for VI as 
well as the LSDA results of Ref. |jl^ yield a much lower 
ratio. 

A comparison with experiment of the calculated spec- 
tra reproduces in each single case some of the experimen- 
tal features but misses others. A plausible explanation 
is that the sample contains a mixture of high and low 
barrier rings. In Fig. 6 we display a superposition of the 
VI and V3 spectra in comparison with the available ex- 
perimental points, assuming that the system takes the 
theoretical ground state at each magnetic field, although 
as already noticed in many cases the triplet-singlet en- 
ergy difference is quite small. With this interpretation, 
the circles are attributed to high barrier rings while the 
triangles and rombuses would correspond to low-barrier 
ones. Although the positions of the experimental peaks 
are not perfectly reproduced (this would involve a fur- 
ther optimization with respect to loq and Rq) the qual- 
itative agreement is good. A richer mixture including 
intermediate-barrier rings would induce a blurring of the 
different branches with the introduction of additional (in- 
termediate) peaks, but would also yield the above discus- 
sion and comparison with experiment more involved. The 
crosses in Fig. 6 are not reproduced by any of our ring 
calculations and, as hinted by the experimental group, 
they could be due to pure quantum dots in the sample. 
We also mention that in the experiment |^ an insufficient 
signal to noise ratio at low energies does not allow the 
detection of the lower branch at w < 10 meV. 

B. Non-interacting spectra 

Figure 7 shows the singlet and triplet spectra when the 
Coulomb interaction is neglected as in subsection II. C. 
Therefore, they correspond to pure particle-hole transi- 
tions in the single-particle level scheme. On a first look 
they seem quite similar to the spectra of Fig. 5. How- 
ever, after a second inspection important differences ap- 
pear. A dramatic one is the softening, i.e., a decrease in 
energy, of the lowest branch in all three cases. In fact 
for _B w 8 T the lowest state has an almost vanishing 
energy. This fact is explained as the formation of quasi 
Landau bands, even for such a small system, with the 
property that single-particle states of the same band are 
quasi-degenerate and therefore, particle-hole transitions 
are nearly gapless. 

There are also sizeable shifts of peak positions and, in 
general, a higher fragmentation is present in the single 
particle picture. This is easy to understand since the col- 
lectivity asociated to the Coulomb interaction is absent 
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in the second case, and therefore the description of collec- 
tive magnetoplasmon states is quite rough. Altogether, 
the correct qualitative result must be understood as an 
indication that interaction effects are not overwhelming 
for the spectrocopic features of a two-electron system. 

IV. MEAN FIELD THEORY AND 
BROKEN-SYMMETRY SOLUTION 

Usually, the exact solution of the many-body 
Schrodinger equation is not possible and one must con- 
tent with approximate solutions. In this sense, the mean 
field aproximation allows in many cases the description 
of relevant physics. Although in principle the mean field 
provides an accurate picture for a big enough number of 
particles, it is interesting to compare it with the exact 
solution in the present case for two particles. This will 
provide a strong test of the mean field result, that will 
allow us to quantify its validity. We remark that for sys- 
tems confined in an external potential such as atomic or 
quantum dot electrons the restriction on having a suffi- 
cient number of particles in order to develop the mean 
field is not as strict as in self-bound systems, such as liq- 
uid drops or atomic nuclei. This is easy to understand 
since in confined systems the fixed external potential is 
already an important contribution to the mean field and, 
therefore it is not solely based on particle-particle inter- 
actions. 

Imposing circular symetry the mean field theory has 
indeed been applied to medium-large quantum dots 
(TV « 20 - 200), in the Hartree [|T||, Hartree-Fock |l|] 
and LSDA schemes [^,^. Nevertheless when using 
symmetry-unrestricted approaches, even at the mean 
field level, it is normally necessary to restrict to smaller 
sizes (N ~ 10) because of the computational require- 
ments H-H. 

A. LSDA ground state 



The overall evolution of the total angular momentum is 
also well reproduced in MF, especially for the triplet. 
Quite strikingly, however, the angular momentum of the 
singlet does not show the step-like evolution peculiar of 
the exact solution. In fact it is taking in most cases clear 
non-integer values. This must be attributed to a sponta- 
neous breaking of the rotational symmetry in the frame 
of the mean field. 

The symmetry of the MF solution is more easily appre- 
ciated in Fig. 9, that shows the MF ground state density 
P = + Pl and magnetization m = P] — Pi, where p^ , 
Pl are the two spin densities, for VI and V3 at _B = 0. 
Both cases correspond to singlet solutions but, while VI 
has a circularly symmetric density and vanishing magne- 
tization, V3 is showing a separation of spin up and down 
densities, leading to an oscillation of the magnetization 
along the ring perimeter (a spin density wave) as well as 
of the total density (a charge density wave). This be- 
haviour is in fact common to all singlet states at B > 
for which (L^) takes non-integer values. This peculiar 
states were predicted by Reimann et al. [26| in narrow 
rings with a bigger number of electrons, and have been 
also studied by us in Ref. pTf . 

The formation of broken-symmetry solutions within 
MF approaches is a well known phenomenon in nuclear 
physics In this field, it has been shown that the MF 
is in fact providing the system structure in an intrinsic 
reference frame, given by proper intrinsic coordinates, 
while there is degeneracy with respect to the remain- 
ing collective coordinates. In a simplifying picture this 
corresponds to assume a certain structure of the system 
relative to its mean field, that in turn may be moving 
thus leading to a restoration of the exact symmetry. In 
fact, as shown in Fig. 1, when performing an angular av- 
erage of the MF density for V3 an excellent agreement 
with the exact result is obtained. As we discuss further 
in the next subsection, the MF symmetry breaking can 
be interpreted as an electronic crystallization, i.e., the 
formation of a Wigner molecule. 



We have applied the LSDA symmetry-unrestricted ap- 
proach developed by us in Ref. |Q to the present two- 
electron rings. Figure 8 summarizes the ground state 
properties, for both singlet and triplet for VI, V2 and 
VS. The comparison with the energies for the exact so- 
lution (Fig. 2) prove that the the mean field is providing 
a quantitatively good description of the ground state en- 
ergy, with relative energy differences around 2%. We 
notice that the MF description of the triplet is slightly 
better than that of the singlet. In fact, after the first 
singlet-triplet crossing the MF fails to reproduce a sin- 
glet state more bound than the triplet one. This can be 
attributed to a more important effect of correlations due 
to the interaction in the singlet than in the triplet which 
are only approximately taken into account within LSDA. 



B. The Wigner molecule 

The symmetry breaking of the MF solution for V3 is 
pointing an incipient electron localization in the intrin- 
sic frame of this system. The formation of these local- 
ized structures has been recently discussed using exact 
solutions by Yannouleas and Landman for two electron 
parabolic dots [TT| and, for rings with 6 electrons by 
Koskinen et al. |12|. These authors have computed the 
exact excitation spectrum and shown that it adjusts to 
rotor bands, thus revealing the existence of the electron 
molecule. 

Even when the exact wave function is known it is not a 
trivial task to notice the existence of a Wigner molecule 
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in the intrinsic frame. It can be hinted by means of intrin- 
sic distribution functions, hke the conditional probabihty 
density p{r2,92i\ri — S,), which gives the probabihty of 
finding an electron at (r2, O21) when the other one is fixed 
at ri = ^. Therefore, it provides valuable insight on the 
intrinsic structure of the two-electron system. In Ref. 
p9[ Berry reviews the use of the conditional probability 
density to describe the formation of electron molecules 
in the doubly excited states of the helium atom. In what 
follows we show that a similar analysis can be used for 
the two-electron rings. 

Figure 10 shows the conditional probability density 
when the first electron is placed at its most probable 
value for the radial coordinate in both VI and V3 con- 
fining potentials. Since this position is taken as the ori- 
gin for the angle 6*21 this implies that the first electron 
is placed in those graphs at (^,0), with ^ = 1.5, 1.8 Og 
for VI and V3, respectively. Notice that the two dis- 
tributions yield a peak for the second electron that lies 
opposite to the first one, i.e., at (—^,0). On a first look 
both functions look rather similar although a closer in- 
spection reveals that the inner hole is more pronounced 
in the case of V3. This is in agreement with a stronger 
localization of the second electron in this case. In or- 
der to manifest more clearly this result it is helpful to 
construct the density in the intrinsic frame and com- 
pare it with the MF one. To this end we assume that 
once the most probable relative positions are determined, 
(x, y) = (^, 0) and (— ^, 0), the intrinsic frame density can 
approximately be obtained by centering in each position 
the conditional probability density. This procedure leads 
to the two lower plots of Fig. 10. Quite interestingly, 
the V3 result is clearly breaking the circular symmetry 
and both are in excellent agreement with the mean field 
densities of Fig. 9. 

Although qualitative the preceding discussion supports 
the interpretation that the MF describes the intrin- 
sic structure of the ring. Consequently, the existence 
of a high barrier in V3 produces an incipient Wigner 
molecule. However, it remains as an interesting task for 
the future to obtain the full set of excited states and dis- 
cuss its rotor-like character, following Refs. jll|,|l^ as well 
as the analytical expressions of Ref. H . 



C. Unrestricted time-dependent LSDA (TDLSDA) 

Based on the mean-field picture, the TDLSDA al- 
lows to describe the collective dipole oscillations. These 
scheme has been used to describe both spin and density 
excitations in quantum dots and rings in a symmetry- 
restricted approach |3^Jl^ . In Ref. 1 23 1 we described the 
application of TDLSDA relaxing the constraint of circu- 
lar symmetry and using a real-time approach similar to 
that of Sec. III. Within TDLSDA this implies the solu- 
tion of the time-dependent Kohn-Sham equations, which 



are the analog of Eq. ^ within the LSDA. We have applied 
this technique to the present two-electron rings and show 
the corresponding singlet and triplet spectra for VI, V2 
and V3 in Fig. 11. 

The TDLSDA provides a rather good description of 
the rings' dipole spectra, as compared to the exact ones. 
We notice that the overall evolution with B and the lo- 
cation of the strength is quite well reproduced for each 
branch. However, finer details of the spectra, such as 
the fragmentation structures are not reproduced in many 
cases. Although based on a single-particle picture the un- 
restricted TDLSDA corrects the lower branch instability 
found in Fig. 7. This is a manifestation of the interac- 
tion effects taken into account by the TDLSDA. From 
the present results we conclude that, despite the present 
rings only contain two electrons, the TDLSDA provides 
a rather sensible picture of the dipole oscillations. 

V. CONCLUSIONS 

In this work we have theoretically analyzed two- 
electron rings motivated by recent measurements of their 
FIR absorption. We have found that in order to explain 
the experimental spectra it is necessary to assume that 
the sample contained a mixture of rings with different 
barrier heigths. As expected, high barrier rings are char- 
acterized by a clear violation of Kohn's theorem, showing 
a sizeable amount of strength at energies uj > 15 meV for 
B < 6 T. On the other hand, low barrier rings approxi- 
mately fulfill Kohn's theorem, with an effective ujq, and 
have lower energy excitation branches more similar to 
the measured ones. A mixture of these two types of rings 
has been used to explain the main features of the exper- 
iment. A more detailed analysis would require includ- 
ing intermediate barrier rings (ideally knowing the actual 
barrier- height distribution in the sample), for which the 
spectra lie in between the preceding two extreme situ- 
ations. Despite the very small number of constituents, 
these systems possess a rich FIR spectrum with an in- 
trincate i3-dependence. 

The evolution of angular momentum and total spin 
with barrier height and magnetic field has been discussed 
and contrasted with independent particle predictions. In 
agreement with previous studies we find that interac- 
tion effects become less important for increasingly narrow 
rings, obviously reflecting a bigger electronic separation 
than in a compact structure and a corresponding increase 
in angular momentum. 

The comparison with the symmetry-unrestricted 
LSDA has allowed us to discuss the formation of an incip- 
ient Wigner molecule for the V3 case and, in general, has 
shown the good quality of this approximation. We have 
proposed a method to obtain the density in the intrin- 
sic reference frame, based on the conditional probability 
density, that quantitatively agrees with the LSDA result. 



6 



The FIR absorption in TDLSDA provides a good overall 
strength distribution and dispersion with magnetic field, 
although it misses some finer details of the fragmentation 
patterns. 
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FIG. 1. Electron density for the three rings (symbols) as a 
function of r. Notice that the exact solution always yields cir- 
cularly symmetric densities. In order of decreasing density at 
r = the results correspond to VI, V2 and V3, respectively. 
Also shown in this plot are the LSDA densities (solid lines), 
circularly averaged in the case of broken symmetry (V3). In- 
set: External potential V{r) for rings VI (solid), V2 (dash) 
and V3 (dot-dash). 
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FIG. 2. Evolution with magnetic field of the lowest singlet 
and triplet states for the three rings. Left panels show the 
state energy (in H*) and right ones its corresponding total 
angular momentum in the perpendicular direction (Lz), in 
units of h. 




FIG. 3. Dependence of the single particle energy levels with 
B for the three potentials. The different branches, taken in 
vertical order at B = IT, correspond for all three cases to 
{n£) = (10), (1 -1), (1 1), (1 -2), (1 -3), where n (radial) and 
£ (^-component of angular momentum) are the usual quan- 
tum numbers. Each branch presents a small spin splitting for 
increasing B, with the spin up states placed below the spin 
down ones. 
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FIG. 4. Same as Fig. 2 for the non-interacting case. 



FIG. 5. Evolution of the FIR spectrum for singlet and 
triplet ground states in the three rings. The results for VI, 
V2 and V3 are displayed in the three rows from top to bottom 
while left and right columns correspond to singlet and triplet 
states, respectively. For comparison, the dashed lines show 
the analytical predicition from Kohn's theorem (j^, see Sec. 
III.A. 



FIG. 6. Superposition of VI (solid lines) and V3 (dashed) 
spectra at different magnetic fields. The experimental data 
are taken from Lorke et al. [9]. See text (Sec. III.A) for 
details. 



FIG. 7. Same as Fig. 5 for the non-interacting case. 
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FIG. 8. Same as Fig. 2 within the LSDA. 



FIG. 9. Results obtained within the LSDA for S = T: 

Lower plots display the density p and magnetization m for V3 
while the upper one shows the density for VI. In the latter 
case the magnetization vanishes everywhere. 



FIG. 10. Upper plots: conditional probability densities 
(CPD) p(r2,^2i|ri = ^) obtained from the exact wave func- 
tion at S = T. The first electron has been fixed at (^,0), 
indicated by a solid dot, where ^ = 1.5, 1.8 is the most 
probable radius for a single electron, in VI and V3 respec- 
tively. Lower plots: Intrinsic-frame density built from the 
upper conditional probability densities as explained in the 
text (Sec. IV.B). 



FIG. 11. Same as Fig. 5 within the LSDA. 
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